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Abstract 

Using numerical simulations, the stability and scattering proper¬ 
ties of the 0(3) model on a two-dimensional torus are studied. Its 
solitons are found to be unstable but can be stabilized by the addition 
of a Skyrme-like term to the Lagrangian. Scattering at right angles 
with respect to the initial direction of motion is observed in all cases 
studied. The model has no solutions of degree one, so when a field 
configuration that resembles a soliton is considered, it shrinks to be¬ 
come infinitely thin. A comparison of these results with those of the 
model defined on the sphere is made. 


1 Introduction 


The non-linear 0(3) model in two dimensions appears as a low di¬ 
mensional analogue of non-abelian gauge field theories in four dimen¬ 
sions. This analogy relies on common properties like conformal in¬ 
variance, non-trivial topology, existence of solitons, hidden symme¬ 
tries and asymptotic freedom. Among various applications, the 0(3) 
model has been used in the study of the quantum Hall effect and in 
models of high-T c superconductivity. In solid state physics it arises 
as the continuum limit of an isotropic ferromagnet, and in differential 
geometry the solutions of the model are known as harmonic maps. 

The classical 0(3) model defined on the sphere, or on compactified 
plane, has been amply discussed in the literature §0- Its finite- 
energy solutions are static solitons or instantons describing lumps of 
energy. The model in three dimensional space-time is not integrable, 
and so to study the time evolution of its solitons one must resort to 
numerical simulations. The stability of the solitons has been analyzed 
in reference fj . Due to the conformal invariance of the theory in two 
spatial dimensions, these soliton-like lumps are unstable in the sense 
that they change their size under any small perturbation, either ex¬ 
plicit or introduced by the discretization procedure. When the lumps 
collide head-on, they scatter in general at 90° to the initial direction of 
motion in the center-of-mass frame 11 - Some time after the collison 
the instability of the solitons makes them shrink at an ever-increasing 
rate and they become so spiky that the numerical simulations break 
down. However, it has been shown i 1 that this instability of the 
model can be cured by the addition of two extra terms to the Lagra- 
gian. The first one resembles the term introduced by Skyrrne in his 
nuclear model in four dimensional space-time [0], and the second one 
is an additional potential term. The solitonic solutions of this mod¬ 
ified model are stable lumps (skyrmions) which repel each other and 
scatter at 90° when sent towards each other with sufficient speed. 

In the present paper we study the evolution properties of the 0(3) 
model when periodic boundary conditions are imposed; this amounts 
to defining the classical model on a two-dimensional torus. This ap¬ 
proach looks more physical than the one on the sphere in the sense 
that the solitons are located in a finite volume from the outset. In 
any case, a comparison between both the toroidal and the spherical 
approaches is certainly of interest, if only to check the consistency of 
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the two results. 

In the next section we present the 0(3) model on the torus, and 
explain the numerical set up in the following one. Solitons of degree 
one in both the 0(3) scheme and its Skyrrne version are discussed in 
section [|, whereas their scattering is studied in section |5| Section ||] 
completes our paper with some conclusions. 


2 The 0 ( 3 ) model on the torus 

The non-linear 0(3) model involves three real scalar fields ^>(x^) = 
{c t> a (x A ‘), a = 1,2, 3} with the constraint that for all x^ = (x°, x 1 , x 2 ) = 
(t,x,y) the fields lie on the unit sphere S 2 : 

= 1- (1) 


Subject to the above constraint the Lagrangian density and the cor¬ 
responding equations of motion read 

£ = ( 2 ) 

d^d^cj) + = 0. (3) 

For any value of t, the fields 4> are mappings from the torus T 2 to the 
sphere S 2 , i-e., they satisfy the periodic boundary conditions 

<f>(x + mL,y + nL) = (j>(x,y), (4) 

where m, n = 0,1, 2,..., and the period L denotes the size of the square 
torus. 

It is convenient to describe the model in terms of one independent 
complex field W, related to 4> via 

r ,W + W ,-W + W \W\ 2 — 1 

(P ~ ^\W\ 2 + V l |IT| 2 + 1 ’ \W\ 2 + 1 U 


Introducing complex coordinates z = x + iy and z = x — iy on the 
torus and using the handy notation d z W = W z . d z {d z W) = W zz , etc., 
the equation of motion ([!) for the static field configurations becomes 


T „ 2WW Z W Z n 
Wzg ~ \W\ 2 + 1 “ °’ 


( 6 ) 
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whereas the boundary conditions ([I]) take the form 

W(z + mL + inL) = W(z). (7) 


Whence, our static solitons or instantons are elliptic functions that 
may be expressed as || 


w =a n 


(j(z — Qj) 


Li oiz-bjY 


( 8 ) 


with the zeros ( a,j) and poles (bj) subject to the selection rule 

J2 a j = J2 b j- ( 9 ) 

j =i j-\ 

The complex number A is related to the size of the soliton and k is the 
order of the elliptic function W. The Weierstrass a-function is defined 
on X 2 and satisfies the pseudo-periodicity property || 

a(z+mL+inL) = (—f)( m + ri + mn ) exp(y(m— in)[z+-(m+in)L])a(z). 

L 2 

(10) 

For a square torus we have the so-called lemniscatic case, where a 
posseses a Laurent expansion of the form 

OO 

= X^ 4J+1 ’ ( n ) 

3 = 1 

where the real coefficients Cj depend on L. 

Note that the periodicity of the torus means that our study can 
be simplified to the consideration of the system in a fundamental cell 
delimited by the vertices 


(0,0), (L,0), (L,L), (0,L). (12) 


We are interested in static finite-energy solutions, which in the 
language of differential geometry are harmonic maps T 2 —> S 2 - These 
maps have been extensively studied in differential geometry & 0- 
They are partitioned into homotopy sectors parametrized by an in¬ 
variant integral index Q , the degree of the map, defined as usual by 
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taking a two-form from S 2 to T 2 via the pull back map W*. For a 
given map W : T 2 —» S 2 we can pull back the Kahler form 


lo = 4i 


A d£ 

a+W’ 


teS 2 , 


(13) 


and define 

Q = i / W*u>, ( 14 ) 

F Jt 2 

where the constant C normalizes Q to an integer. Expanding ®) in 
terms of z and setting C = — 87 r we obtain 


Q 


1 f [HT| 2 - |W g | 2 

ttJt 2 (|1F| 2 + 1) 2 


dxdy. 


(15) 


The potential energy V, as derived from the Lagrangian (||), and the 
topological index satisfy 

V > 2ir\Q\. (16) 

The instanton-solutions correspond to the equality in (|l^): Solutions 
carrying Q > 0 (Q < 0) imply W s = 0 (W z = 0), which are the 
Cauchy-Riemann conditions for W being an analytic function of z 
(z). 

The periodic solitons @ have been studied in a variety of con¬ 
texts. In reference [[0]], for example, they were used to compute the 
contribution of instantons to the partition function. 


3 Numerical procedure for the time 
evolution 


So far we have discussed the static field configurations. Now we 
concentrate on their dynamics, paying particular attention to their 
stability, scattering properties, etc.. As our model is not integrable, 
the study of the evolution of our fields requires numerical techniques. 
Hence, we treat configurations (||)-(j9|) as initial conditions for our 
evolution, studied numerically. The field cj) is evolved according to 
the equation of motion ([ 3 ]). We compute the series (IT) up to the 
fifth term, the coefficients Cj being in our case negligibly small for 
j > 6. The numerical set up is similar to that used in our previous pa¬ 
pers (see [12] for example): We employ the fourth-order Runge-Kutta 
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method and approximate the spatial derivatives by finite differences. 
The Laplacian is evaluated using the standard nine-point formula and, 
to further check our results, a 13-point recipe is also utilized. We work 
on a 200x200 periodic lattice ( n x = n y = 200) with spatial and time 
steps 5x=5y= 0.02 and <5f=0.005, respectively. The size of our torus is 
then L = n x x 5x = 4. 

Unavoidable numerical truncation errors introduced at various 
stages of the calculations gradually shift the fields away from the unit 
sphere ([!]). So we rescale 




every few iterations. The error associated with this procedure is of 
the order of the accuracy of our calculations. Each time, just before 
the rescaling operation, we evaluate the quantity y = (f).(f> — 1 at 
each lattice point. Treating the maximum of the absolute value of y 
as a measure of the numerical errors, we find that max|/r| ~ 10~ 8 . 
This magnitude is useful as a guide to determine how reliable a given 
numerical result is. Usage of an unsound numerical procedure like, 
say, taking Sx < 5t in the Runge-Kutta evolution, shows itself as a 
rapid growth of max|/i|; such increase also occurs when the solitons 
become infinitely spiky. 

We use the global U( 1) symmetry of @ to choose A real; the value 
A=1 has been used in all our simulations. 


4 Solitons of degree one 

4.1 0(3) case 

From the theory of elliptic functions we know that the simplest non¬ 
trivial elliptic functions are of order two. This implies that the 0(3) 
model on the torus possesses no single-soliton solutions. This fact 
may also be understood in the context of differential geometry The 
harmonic maps M —» S 2 (M an orientable surface) have holomorphic 
representatives (instantons) of any degree provided that it is greater 
than the genus of M || |H]]. Clearly, for M = T 2 the index of the 
maps must be greater than unity. Note that the degree Q in ( |l5| ) is 

1 We thank J.M. Speight for showing us a simple mathematical proof of this fact. 
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numerically equal to the order k of (|j) only when the latter is greater 
than one. Thus, an order-one solution (the trivial solution) carries 
degree zero, not one. 

In order to study a single soliton on T 2 , we ignore the selection 
rule (3) and take the configuration 


W\ 


<j(z — a) 
<t(z — b ) ’ 


a 7^ b, 


(17) 


which describes a quasi-periodic lump that instead of (|?[) satisfies 

Wi(z + mL + inV) = exp [y (m — in){b — a)]W\(z). (18) 

1j 

But a periodic solution may be constructed by taking a field whose 
values in the sub-cell of vertices 


(1,1), ( L-l,l ), (L-l,L-l), (, l,L-l), 1<&L, 


are given by W\ and in the rest of the fundamental cell (12) are given 
by a suitably chosen interpolating function. 

So let us periodize W\ along the x-axis with the help of the ansatz 


Wh = A(y) tanh[a(x — L)\ + B(y)\ x e [L — l, L + l], y e [0, L\. (19) 


(note that [L — l, L + l\ = [0, l] U [L — l, L\). The complex functions 
A(y) and B(y) are obtained by demanding periodicity and continuity 
of© and ©■ One deduces 


A(y) = \W\(l + iy) — W\(L — l + zy)]/2tanh(ai), 
B(y) = \Wi(l + iy) + W\(L — l + iy)\/2. 

Therefore, our horizontally-periodic configuration is 


W H {x,y) 


Wi; x e[l, L — l], ye[0,L\; 

W h - xe [L-l,L + l], ye [0,L]. 


( 20 ) 


( 21 ) 


A similar periodization is now performed on Wh along the vertical 
axis. It turns out that the field W p thus obtained is periodic in both 
x and y. It has the appearance 


7 





W H ; xe [0, L\, ye[l,L — l]; 

C{x) tanh[a(?/ — L)\ + D(x)\ xe[0,L\, ye[L — l,L + 


l], 

( 22 ) 



with 

C(x) = {Wh(x + il) -Wh[x + i(L - l)]}/2tanh(al), , . 

D(x) = {W H {x + il) + W H [x + i(L-l)}}/2. 1 ’ 

For the length l we may take ten lattice points, so that l = 0.2<L = 
4. We set the value of a equal to 20 , and for the zero and pole of ( [I?] ) 
we elect 

a = (2.05,1.75), b = (1.95, 2.25). (24) 

We have numerically checked that the ansatz (^) has Q = 1, and so 
it may be regarded as a map T 2 —»• S 2 of degree one. 

The upper half of figure (|l|) illustrates the periodization of W\ for 
a representative line of the fundamental cell. The lower half exhibits 
the total energy density associated with the periodized field W p . It is 
apparent from this picture that the periodization procedure introduces 
some perturbation at the borders of the grid in the form of small 
folds. Under the numerical evolution these perturbations propagate 
towards the center of the lattice and collapse the lump. In order to 
minimize such an effect, we improve the initial conditions by ironing 
out the folds. We do this by implementing a damping function, 7, that 
rescales dfp p —> 7 dt4> p , 7 < 1. The absorbtion is switched off at the 
time (to) when the folds have disappeared; the resultant configuration 
serves us as a better, improved set of initial conditions. 

During the preparatory stage the total energy undergoes a small 
decrease, in conformity with the absorption that is taking place. Once 
the latter is turned off, the energy settles near the expected value of 
one and remains constant until the time ( tf ) when the total energy 
density becomes so spiky that the numerical procedure breaks down 
[see figure [|. Moreover, we have checked that these results do not 
depend on how the initial conditions were prepared, nor on whether a 
9-point or a 13-point laplacian operator was used in the simulations. 
Having performed many such simulations we are convinced that our 
results are genuine, i.e., the shrinking is genuine and not a numerical 
artifact. 
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4.2 Skyrme case 


Next we look at possible ways to stabilizing our solitons. Guided by 
the experience with the 0(3) model in the plane (where the shrinking 
can be prevented by the addition of two extra terms to the Lagrangian 
-the Skyrme and the potential term-) we consider the possibility of 
adding the Skyrme term alone. Adding such a term [|li| to the 0(3) 
Lagrangian we have: 


skyrme — ^ ■ (<9^0) 




or, in the simpler CP 1 formulation 
„ IkLil 2 — 2\W Z \ 2 i an 

^skyrme ^ i I T/l/12 \ 2 ' ^1 


(i + l^l 


|Wz|2 {\w t \ 2 -\w z \ 2 ). 


(i + W 


2\A ’ 


(25) 

(26) 


The associated static equation of motion reads 


o = w z ~ - 


2 WW z Wg 

\W\ 2 + 1 


+ (\nr\ if- [2W ^ w ‘ w ’ ~ W ^ W P 
+w gg w g w z - W zg ( |VP0| 2 + \w- z \ 2 ) 


+ 


2 W 


\W\ 2 + 1 


(|w / 2 | 2 -|w= 


2\2i 


W gg (W g ) 2 + w zz w z w s 

(27) 


Indeed, we find that thanks to the extra term the energy density of 
the lump does not increase indefinitely, but instead it oscillates with 
time in a stable manner. In figure (j|) we show the evolution of the 
amplitude of the total energy density for 0\ = 0.001. Qualitatively 
similar pictures are obtained for values of 9\ as small as ~ 0.00015; 
for smaller values W p is no longer stable. 

Note that W p does not exactly satisfy the equation of motion (|27|), 
for the term 


46»i 

(|IL| 2 + 1) 2 


[W zz (W z f 


O W|^| 4 
|IL| 2 + 1 


does not vanish. Nevertheless, the smallness of 6\ means that our 
Skyrme model is only a slight perturbation of 0(3), and hence W p is 
a good, if approximate, solution. 
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Worthy of remark is the fact that (^5|) does not require a potential¬ 
like term to stabilize the lumps. In the 0(3) system in the plane with 
fixed boundary conditions, by contrast, such a term was needed to 
prevent the solitons from expanding. 


5 Solitons of degree 2 


5.1 0(3) case 

We now move on to the interesting question of collisions, limiting 
ourselves to two solitons. It is important to bear in mind that the 
preparatory stage devised for the pathological single-soliton case of 
section [| is not required for lumps of degree > 2. The initial field is 
given by an order-two fuction of the form ®-i): 


w. 


cr(z — ai) a{z — a 2 ) 
a(z - 61 ) a(z - b 2 ) ’ 


ai + a 2 = h + b 2 . 


(28) 


First, consider the situation when the solitons are symmetrically 
positioned along the horizontal axis and boosted towards each other 
with relative velocity v = (0.2, 0). We select the zeros and poles to be: 


01 = (0.77,1.95), o 2 = (3.25,1.95); 
61 = (1.32,1.95), b 2 = (2.70,1.95). 


The solitons gradually shrink and then undergo a gradual expansion 
as they approach each other. They collide at the centre of the grid 
and merge into a complicated ringish structure, where they are no 
longer distinguishable. After this process the solitons get narrower 
and narrower as they re-emerge at right angles to the initial direction 
of motion. Due to their instability, the shrinking process goes on until 
the solitons get so spiky that the numerical procedure is no longer 
reliable; this occurs for t ~ 7, when max|/i| as defined in section [| 
reaches ~ 10 -4 and higher. 

A numerically interesting feature of the periodic 0(3) model is that 
the scattering can also be observed when the solitons are sped ‘away’ 
from each other, towards the borders of the fundamental cell. This is 
a good way to test the correctness of our periodic lattice. Applying 
this to the solitons defined by (|29| ) we also observe the scattering at 
90°. A representation of this process can be viewed in figure (||). 
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A typical head-on collision with the solitons initially placed along 
a diagonal is illustrated in figure (j5|). The initial position is achieved 
by the arrangement 


a\ = (0.95,0.75), a 2 = (3.05,3.25); 
b x = (1.22,1.95), b 2 = (2.78,2.05). 


(30) 


After boosting the solitons away from the centre with initial veloc- 


|u|=0.2-, the lumps collide at the corner (0,0)=(4,4) 


ity v = t?(M) 

and re-appear from (0,4) = (4,0) at right angles to the initial direction 
of motion. Of course, all four corners are nothing but the same point; 
there the lumps meet, coalesce and scatter off as already explained. 
Shortly afterwards, the instability of the system manifests itself in the 
usual manner, as reflected by the 0(3) curve in the graph Emax{t) 
in figure (5j). This diagram also includes the resulting curve of the 
Skyrme version, as described in subsection ( |5.2p below. Also, when 
situated in an arbitrary, non-symmetrical way within the fundamental 
lattice, the solitons always scatter at ninety degrees when sent head- 
on against one another (we discuss in more detail this situation in the 
next section). 

We may interpret the instability of (28) under numerical simula¬ 
tions as follows: The solitons start off satisfying the selection rule 
a\ + a 2 = 61 + 62 , which links them in some manner. Due to inevitable 
round-off errors during the numerical simulation, the field gets per¬ 
turbed and so is only approximately described by the original field 
configuration. As the perturbation is quite small it will excite mainly 
the degrees of freedom which are zero modes of the original configu¬ 
ration. Thus, in particular, aj and bj will start evolving but in order 
to remain close to the original configuration they will keep the con¬ 
straint unbroken. Such evolution may lead to aj and bj, pairwise, 
coming close together. This corresponds to the solitons shrinking. To 
see this note that |aj — bj\/2 determines the size of the j’-th soliton. 
Note that this shrinking is essentially of the same type as the well 
known shrinking of any number of solitons on the sphere. We would 
like to stress that since analytical solutions exist in all topological sec¬ 
tors of index > 2 , this lack of stability of our two-soliton system is of a 
different nature than the instability of the single-soliton configuration 
(and so non-existence of a one-soliton static solution) discussed in the 
previous section. There the solution does not exist on the lattice or in 
the continuum; here the solutions do exist in the continuum but are 
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unstable and putting them on the lattice introduces a perturbation 
which sets off the instability. 


5.2 Skyrme case 


Let us now consider the Skyrme Lagrangian (26) as applied to two 
solitons. Head-on collisions along the horizontal axis corresponding to 
the set up ( |29|) proceed as in the pure 0(3) scheme. The Skyrme term, 
however, prevents the lumps from shrinking indefinitely and renders 
them stable; their motion can now be followed for as long as desired. 
For instance, the skyrmions proceed as in figure (|4|) but, after 90° 
scattering at the lattice point (0,2)=(4,2), they continue their journey 
and collide thrice more, reach again their t =0 positions and proceed to 
repeat this cycle anew, as suggested by figure ®. Note the coalescence 
of the lumps in the corners (subplot t = 11) which, as mentioned 
before, are nothing but one and the same point. The plot of the 
corresponding Emax(t) -not shown- is very much like the one drawn 
in figure (^) [dashed curve]. All two-skyrmion cases shown in this 
paper correspond to 6\ = 1/2000, but the same qualitative behaviour 
is found for values down to ~ 0.00007. Smaller values cannot prevent 
the lumps from getting too thin, leading to the breakdown of our code. 

An example of solitons located at two arbitrary cell points is given 
by the parameters 


ai = (0.77,1.30), 
b\ = (1.32,1.95), 


a 2 = (3.25,2.70); 
b 2 = (2.70,2.05). 


(31) 


When these solitons are sent to collide head-on, with or without a 6\ 
term, the scattering, as usual, takes place at 7 t/2 radians. We shall 
depict this event within the stable format of the modified model. 

Figure ( 0 ) refers to two skyrmions directed towards each other with 
v = (0.18047,0.0862), the speed being 0.2. The coordinates ( x,y ) 
correspond to the position of the amplitude Emax{t). The labels A-E 
are a guide as to the path followed by one of the lumps, the route 
of the other being given by the corresponding symmetrical points. 
So a skyrmion-lump starts at A and after 90° scattering around the 
centre it continues its itinerary to the position B , where it disappears 
to re-emerge at C. Thence the extended structure heads south-east 
and, having reached point D at to « 14.5, it suddenly changes its 
path to move south-west (point E), unequivocally signalling that a 
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second 90° scattering has taken place. Regarding the other colliding 
entity, the one starting at « (1.5,0.9), we can see that the said second 
collision changes the trajectory of this skyrmion from the north-west 
to the north-east direction. Our numerical simulation terminates at 
tp « 30, where F denotes the end of the leg started at E. 

Note that it is also possible to imagine our solitons as evolving on 
the surface of a doughnut in 71$, obtained by rotating the circle of 
radius r and circumference L (the size of the flat manifold To) about a 
coplanar line [Z axis, say) that does not intersect it. The coordinates 
(x = T, y = d) serve as the angle of rotation of the plane of the 
circle and the angle on the circle itself, respectively. The parametric 
equations of such a torus are the standard 

X = [R + r 003 ( 7 ?)] cos(T), Y = [R + r 003 ( 7 ?)] sin(T), Z = r sin($). 

(32) 

Both the radius r and the distance R from the centre of the circle to 
the axis of revolution (Z) can be calculated from L. With the help 
of (|32|), the distance d of a lump on the surface of the torus from the 
origin (0,0,0) may be computed via d = \f X' 2 + Y 2 + Z 2 . 

Finally, regarding the situation where the initial velocity of the 
lumps equals zero, we recall that in the Skyrme model on the sphere 
the solitons (represented by an approximate held solution) slighty 
move away from each other, showing the presence of a repulsive force 
between them 0- On the other hand, on the torus we have found 
that our skyrmions, also an approximate solution of the equations of 
motion, undergo no translation at all as the time elapses. This agrees 
with our expectations: The net repulsive force on a given lump is now 
zero due to the presence of similar entities in neighbouring lattices. 


6 Concluding remarks 

With the help of numerical simulations we have investigated some 
stability and scattering properties of the non-linear 0(3) and Skyrme 
models with periodic boundary conditions in (2+1) dimensions. 

The toroidal 0(3) theory has the distinctive feature of possess¬ 
ing analytic soliton solutions only of degree two and higher. This is 
because the defining fields are elliptic functions. We studied a single- 
soliton case through a periodic ansatz, which has turned out to be 
unstable: After some time the lump of energy grows too spiky and 
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the numerical procedure breaks down. Since there are no analytical 
solutions with topological charge equal to unity, we may regard the 
above instability as intrinsinc to the model rather than an artifact 
of our numerical method. However, our ansatz has become stable 
upon the addition of an extra term (a Skyrme-like term) to the 0(3) 
Lagrangian and, remarkably, under such circumstances our proposed 
field serves as a good, if approximate, soliton solution of degree one on 
the torus. Note that unlike the more familiar model on compactified 
plane, where a second extra term is required to stabilize the solitons, 
our model necessitates only a Skyrme term to achieve so. In this 
sense, our model resembles more closely its parent (3+1) dimensional 
version, where no second extra term is needed, either. 

With regards to collisions we have limited ourselves to those in¬ 
volving two solitons (analytic solutions of which do exist on T 2 ). Scat¬ 
tering at 90° was observed in all cases considered, both in the pure 
0(3) and Skyrme schemes. Within the framework of the former the 
lumps shrink unstably as they evolve, behaviour set off by the per¬ 
turbation brought about by the discretization procedure. Such in¬ 
stability is of the same type as that exhibited in the familiar 0(3) 
model defined on the compactified plane or topological sphere. As in 
the one-soliton case, the sole addition of a Skyrme term stabilizes our 
two-soliton system. Finally, when the skyrmions start off from rest 
they remain motionless under the numerical simulation, in contrast 
with the skyrmions on the sphere where they move away from each 
other due to a repulsive force amongst them. 
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Figure 1: Above : W\ (dashed line) and its periodized version W p (solid line) 
along the line y — 2 of the fundamental cell. Below : A typical picture for 
the total energy density associated with W p at t = 0. The folds at the edges 
are caused by the periodization procedure. 




Emax 





Figure 2: Above : Total energy density at t 0 —' 0.8, corresponding to our pre¬ 
pared. improved initial one soliton-like configuration. Below : The maximum 
value of the total energy density (Emax) and the total energy vs. t. The 
lump grows infinitely tall soon after tf ~ 3.5, and the numerical procedure 
collapses. 
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Emax=118.2 



Emax=83.37 

t=3.5 



Figure 4: Total energy density corresponding to 0(3) solitons moving away 
from the centre with v = (0.2,0). The amplitude of the lumps gradually 
decreases as they approach each other, reaching a maximum when they coa¬ 
lesce. 
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Emax=148.8 



Emax=553.7 



Figure [I] (continued): The solitons scatter at 90°. They become very spiky 
as time progresses but, as shown in figure (||), this is corrected by adding a 
Skyrme term to the Lagrangian. 


20 







Before (B) 


After (A) 





Figure 5: 0(3) solitons moving away from the centre along the (0,0)-(4,4) 
diagonal (B). They collide at the corners and scatter at right angles (A). 
When the model is supplemented by a Skyrme term the lumps are stable, as 
shown in the accompanying graph Emax{t). 
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Emax=23.92 ; t=11 Emax=43.13 ; t=17.5 



Figure 6: Skyrmion scattering {9\— 1/2000) for v = (0.2,0). After scattering 
like the 0(3) solitons of figure 4, the skyrmions do not collapse but go on 
to collide at £=11,17.5,25.5 and so forth. At every occasion they scatter at 
right angles. This cycle repeats itself indefinitely. 
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Figure 7: Trajectories of the position of Emax corresponding to the head-on 
scattering for skyrmions arbitrarily situated in the basic cell. The labels A-E 
indicate the itinerary of one of the lumps. 
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